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ABSTRACT 

Spiral galaxies are observed to exhibit a range of morphologies, in particular in the 
shape of spiral arms. A key diagnostic parameter is the pitch angle, which describes 
how tightly wound the spiral arms are. Observationally and analytically, a correlation 
between pitch angle and galactic shear rate has been detected. For the first time, we 
perform a suite of N-body simulations to calculate and compare the pitch angles of 
both individual density waves and overall spiral structure by use of two independent 
techniques. We find that higher galactic shear rates produce more tightly wound spiral 
arms, both in individual mode patterns (density waves) and in the overall density en- 
hancement. Although the mode pattern pitch angles by construction remain constant 
with time, the overall logarithmic spiral arm winds over time, which is consistent with 
both the observational scatter in pitch angle versus shear seen from observations, and 
the recent idea that multiple mode patterns may interfere with each other to create 
apparently winding spiral arm structures. 

Key words: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 
- galaxies: structure 



1 INTRODUCTION 

Spiral arms are the defining structure of any spiral galaxy. 
Simulations ha v e sho wn that their transient nature (e.g. 
ISellwoodl |201(A 1201 ll ) invokes significant effects on the 
dynamics of baryonic ma tter in the galaxy, parti c ularly 



clouds (jD'Onghia et al.1 l2012h . Reported physical mech- 
anism^) responsible for driving the spiral structure in - 
clude swing ampli fied (| Julian fe Toomrelll966l:lToomrelll98lh 
density waves dLin fe Shul 11964 iGoldreich fe Lvnden-BeU l 



dynamics oi baryonic matter m tne galaxy, particularly 
in the disc plane (e.g . Goldreich & Lynden-Bell 1965; 



Sellwood fe Binnevl 120021 : [Roskar et al. 2008; iGrand et all 



2012a : b), namely the well known radial migration phe- 
nomenon ,_as_wejl_as_Ji^ay^ and possibly disc thicken- 
~vf 120091: lLoebman etaD l201ll : 



ing (e.g. Schonrich & Binne 



19651 : iLvnden-Bell fe Kalnaisl Il972l : Sellwood & Binned 
2002h , multiple spiral/bar wave modes ([Sparke &; Sellwooc 
19871: iMasset fe Tagger] Il997l: IMinchev fe Quillenl 120061 
Minchev fe Famaevll201ol : lOuillen et al.ll201ll: iRoskar et al 
and co-rot ati ng structures (|Wada et al.l l201ll 



ISolwav et al.ll2012l : IMinchev et al.ll2012h . Because all of this 
is directly dependent on the spiral structure, it can be ex- 
pected that different spiral morphologies imprint the dy- 
namical history of a galaxy in different ways and degrees. 

There are several reported causes of spiral struc- 
ture, such as external perturbations from tidal inter- 
wit h 



Grand et al.1 l2012al lbl: iFuiii fe Babal l2012h . Recent obser- 
vations have highlighted the possible existence of both 
densi ty wave and co- rotat i ng type spiral arms in nature 
e.g. iMerrifield et al.1 12005L 120061: iMeidt et al.1 12008L [2 009 ; 



Egusa et al.l 120091: IFovle et al.ll201ll : ISpeights fe Westpfah] 
201 ll : IFerreras et al.l bOli T some of whom find evidence 
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actions with a satellite galaxy (e.g. 
Toomre & Toomrelll972l: [P'Onghia et al.ll2010l : iDobbs et"aD 
20101 : IPurcell et al.l l201lh . The presence of a central bar 
may also be a driver o f density waves and lead to grand 
design spiral galax ies (|Buta et al.l 120051 : ISalo et all 120101 : 
lAthanassoulal l2012l and references therein). Spiral struc- 
ture may also grow from internal disc instabilities such 
as local Toomre instabilities, or arise from the pres- 
ence of a large mass perturber such as giant molecular 
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against the classic spiral density wave theo ry either by 
use of the offset method (jEgusa et al.l 120091 ) that mea- 
sures the offset of s t ar for ming tracers (e.g. IFovle et"aD 
120111 : IFerreras et all [2012) or the some form of the 
Tremaine- Weinberg Equations (e .g. IMerrifield et al.1 12005. 
120061 : ISpeights fe Westpfahj|2uTlh . and some who find ev- 
idence for multiple density waves jMeidt et al. 2008, 2009). 
Whether the mechanism for driving spirals be some form of 
density waves or otherwise, it appears there is some mech- 
anism that shapes galaxy morphology into a preferentially 
tight or open structure. 

A correlation between the shape of the galactic ro- 
tation curve and the orientation (pitch angle) of a spiral 
arm-like density enhancement was derived analytically by 



2 Grand et al 



lJulian fe Toomrd (|l966h (see also iFuchs et al] l2Q05h - who 
used the Hill approximations applied to the shearing sheet, 
in which the evolution of the response density to a large 
perturbing mass is calculated as swing amplification am- 
plifies the density around the perturber. They showed that 
the density enhancement in this context is predicted to show 
smaller pitch angles (i.e. more tightly wound structure) with 
increasing amount of shear present. 

This co r relatio n is supported observationally as well. 
ISeigar et al.l (|20Q5h looked at a sample of several barred 
galaxies and used a Fourier analysis to draw the spiral shape, 
and found evidence for the shear rate dependency of the spi- 
ral arm pitch angle. This s urvey was later ext ended and the 
conclusion strengthened bv lSeigar et al.l (|2006h . 

In this paper, we aim to test both theoretical and ob- 
servational predictions by running a suite of N-body sim- 
ulations, each with varying shear rates. For the first time 
we investigate the pitch angles of individual spiral patterns 
in N-body simulations by isolating the patterns from the 
system and calculating the spiral phase of the ra-th mode. 
Because we assume a density wave nature in order to study 
individual patterns, the pitch angle is assumed to be con- 
stant. We find that the discs of higher shear rate exhibit 
systematically smaller pitch angles as predicted from the 
analytical studies mentioned above. We also trace the over- 
all spiral feature and its pitch angle as a function of time. 
When measured in this way, we find that the pitch angle de- 
creases with time because the appare nt spiral arm feature s 
behave as short-lived transient spirals (Sellwood 2010, 2011), 
that appear to wind in si mulations (|Grand et al.l l2Q12al lbh . 
The observational work of Sei gar et al.l (|2006r h which mea- 
sures the pitch angle and shear rate of many spiral galaxies, 
reveals several different observed pitch angles for a given 
shear rate. We find that this scatter in pitch angle versus 
shear appears to be consistent with the range of pitch an- 
gles that the winding arms in our simulations exhibit. The 
pitch angle evolution might also shed light on whether or 
not the spiral arm can be an apparently winding structure 
owing to the interference of multiple wave mode patterns 
(|Comparetta fe Quillenll2012h . 

The simulations are described in Sec. 2, the analysis and 
results are described in Sec. 3 followed by the discussion in 
Sec. 4. 



2 SIMULATIONS 

The simulations in this p aper are performed with a hierar- 
chical Tree N-body code (jKawata & Gibson! [2003). We per- 
form three simulations, each of which consists of a spherical 
static dark matter halo (and a spherical static stellar bulge 
component in two of the simulations) and a live stellar disc. 

The dark matt er halo density profile follows that of 
INavarro et "all <|l997h : 
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Figure 1. The circular velocity at t = for simulation R (solid), 
F (dashed) and K (dot-dashed) lines respectively. 
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Figure 2. Galactic shear rate for simulation R (solid red), F 
(dashed green) and K (dot-dashed blue) lines respectively. 

ter sphere is equal to 200p C rit (where pcrit — 3Hq /8ttG; the 
critical density for closure): 

r 200 = 1.63 x 10- 2 h^lwc. (2) 

We assume Q 0.266, Q b 0.0044 and H 71 km s" 1 
Mpc -1 . 

The spherical static stellar bulge com ponent is modelled 
by the Hernquist profile (|Hernquistll9 90) . which is described 
by: 



Pcim — 



3Hq Q,q — Q,b 



8ttG tt cx(l + ex) 



(1) 



where p c is the characteristic density described by 
INavarro et al.l (|l997l ). the concentration parameter, c = 
?"2oo /r a , and x = r/r 200- The scale length is r s , and r2oo 
is the radius inside which the mean density of the dark mat- 



( v M b a 



2tt r (r + a) 



(3) 



where Mb is the total bulge mass and a is the scale length. 
The scale length is set to the effective radius, R e = 1.8153a. 
We apply a compacting factor, b, to s cale from the emp irical 
relation of the bulge effective radius (jShen et al J 12002ft : 
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Simulation 


M d (M ) 


^d(kpc) 


M vir (M Q ) 


c 


M b (M ) 


b 


R 


1.0 x 10 10 


3.5 


2.5 x 10 12 


5 






F 


5.0 x 10 10 


3.5 


1.5 x 10 12 


15 


4.0 x 10 10 


0.5 


K 


5.0 x 10 10 


3.5 


1.0 x 10 11 


15 


1.0 x 10 11 


0.01 



Table 1. Table of initial model parameters for simulation R, simulation F and simulation K. 
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Figure 3. Toomre stability parameter, Q, at t = for simulation 
R (solid red), F (dashed green) and K (dot-dashed blue) lines 
respectively. 



Hence, the resultant scale-length is defined by a = 
&i?e ; o/1.8153. 

The stellar disc is assumed to follow an exponential sur- 
face density profile: 

pd >* = a Md p — sech2 ( ~ ) exp \ ~~Er~) • ( 5 ) 

We use N — 1 x 10 6 particles for each simulation to describe 
the disc. We apply a fixed softening leng th for star particles 
with t he spline softening suggested by IPrice &; Monaghanl 
(2007). The softening length for galaxies F and K is 340 pc, 
and for simulation R it is 200 pc. The model parameters 
for the simulations are summarised in Tabled and the ro- 
tation curves are shown in Fig. [T] The rotation curve that 
rises with increasing radius (simulation R) is produced by 
setting a large halo mass with a low concentration param- 
eter, c, in order to extend mass to the outer regions of the 
disc. Because of such a low concentration of dark matter 
mass in the central region, the disc ma ss must be lowered 
in or der to prevent a bar from forming ([Ostriker &; Peebles! 
1 19731 ). In this way, we avoid the added complication of the 
bar component and restrict the study to spiral galaxies only. 
To produce the flat (simulation F) and keplerian-like (sim- 
ulation K) rotation curves, a bulge component is included. 
For simulation K, this is a very compact and massive bulge. 
Although this case is unrealistic, we include it in order to 
emphasise the effect of galactic shear on spiral morphology. 



The radial profile of the galactic shear rate at t — 0, given 
by T = 1 — (R/Vc)(dV c /dR)j for each simulation is shown in 

Fig.!! .. . 

lJulian &; Toomrd (| 19661 ) note that the response density 
can be weakly dependent on the Q stability parameter (de- 
fined as Q = aRhv/3.36GY,). Therefore, we set the initial Q 
value for all our simulations (Fig. [3} to approximately 1 over 
the radial range 4 < R < 12 kpc, so as to minimise any effect 
of this global parameter. 



3 RESULTS AND DISCUSSION 

Here we present the analysis and results of our two tech- 
niques for measuring pitch angles: mode pattern analysis 
and apparent spiral arm peak trace method. An important 
difference between these techniques worth noting is that the 
mode pattern analysis focuses only on individual density 
waves that may be present in the stellar disc, whereas the 
apparent spiral arm peak trace method analyses the spiral 
density enhancement as a whole. 

Before we describe these two analysis techniques, we 
outline the method of pitch angle calculation. First, as- 
sume the positional information of a density enhancement 
is known. Then, the pitch angle can be expressed in terms 
of the position in polar coordinates: 

= arcta n M, (6) 

where the distance, dg, is the spatial distance in the az- 
imuthal direction denned as de — RAO. We assume loga- 
rithmic spiral arms, which means that we can use equation 
((6]) to fit a linear function described by R = bRO + a, where 
b is the gradient of this line, equal to the tangent of the 
pitch angle. The next step is to recover the positional infor- 
mation of individual mode patterns and the apparent spiral 
arm peak in order to apply the linear fitting. 

3.1 Pitch angles of mode patterns 

In this section, we focus on the theory of swing amplified 
modes and their pitch angles. In the context of spiral density 
waves, swing amplification is the mechanism for the growth 
of density centred around a single co-rotation point. Studies 
of the shearing sheet approximation tell us that the orien- 
tation of the response density should depend on the rate of 
shear in the sheet. To extend this to galaxies, it is the mor- 
phology of the density wave that depends on the degree of 
galactic shear as a function of radius. The shape of a wave 
mode pattern is time independent, because each pattern has 
a constant pattern speed, therefore the pitch angle remains 
constant over time. 

In this analysis, we focus on strong patterns because 
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Time(Gyr) Time(Gyr) Time{Gijr) 



Figure 4. The amplitudes of spiral modes m = 1 (red), 2 (green), 3(blue), 4(cyan), 5 (yellow), 6 (black) and 7 (magenta), normalised 
to the axisymmetric m = mode, as a function of time for simulation R (left), F (middle) and K (right) simulations. Vertical dashed 
lines represent the time window of the Fourier transform applied in Section 3.1. Strong contrasts are observed owing to the low Q ~ 1 
parameter that induces spiral structure. 
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Figure 5. Top panels: Power spectra of simulation R for the m = 5 mode (left), F for the m = 3 mode (middle) and K for the m = 2 
mode (right). Prominent ridges (dark pixels) span between 5-10 kpc in most cases. Bottom panels: Dominant mode pattern phase 
positions in the azimuth-radius plane for the corresponding spectra above. Each colour represents a different pattern. The symbol size 
reflects the relative amplitude of the pattern at a given position. The black lines show the lines of best fit for each pattern. 



they give the greatest signal, whereas weaker patterns may 
be confused with noise. In order to find patterns of signifi- 
cant amplitude, we first search for dominant modes i.e. wave 
modes of m spiral arms that exhibit large amplitudes. The 
amplitude of a given wave mode, m, is calculated from the 
quantities: 

N 

W™(r,t) = ^cos(m^), 

i 

N 

WT{r,t) = Y J sHmOi), (7) 

i 

where 9% is the azimuthal angle between the radial vector of 



the particle and a common reference vector. The amplitude 
is then calculated as: 

A m (r,t) = (W?(r,tf + Wr(r,t) 2 ) 1/2 - (8) 

The maximum amplitude in a radial range 4—10 kpc is found 
for modes m — 1 — 7 over the entire 2 Gyr of the evolution 
for each simulation, shown in Fig. U It is evident that there 
are several modes present in each simulation. We focus our 
analysis on the most prominent modes in each simulation, 
i.e. m — 5, 3 and 2 for simulations R, F and K respectively. 

The first step to recovering the positional information 
of a pattern of given multiplicity, m. from the simulation, 
is to compute the power spectra of the strong modes. We 
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compute the power spectra by taking the Fourier transform 
of each component in equation (0 : 

W™(r,u>)= Wr(r,t)e iut h(t)dt, 

WT(r,u>) = ^ zWr(r,t)e^h(t)dt, (9) 

where h(t) denotes the Hanning function. T\ and T2 denote 
the beginning and end of the time window of the Fourier 
transform, which is chosen to be at around a relatively late 
epoch of the simulation when the system is more stable than 
in earlier times (see dashed vertical lines in Fig. [4}. It is also 
centred around a peak of the most dominant mode present 
in each case, and spans At = 256 Myr. 

The amplitude in each frequency as a function of radius 
is then calculated via: 

A m (r,u) = {W?{r,w) 2 + Wr(r, U ) 2 ) 1/2 . (10) 

For each of our galaxies, the power spectrum of the domi- 
nant mode is calculated from the square of the amplitude 
given in equation (|10p , and shown as a function of radius 
and the pattern speed, Q = cj/ra, in the top row of Fig. 
[5] Because simulations generally possess several patterns for 
a given mode that can overlap in radius (e.g. see Fig. 4 of 
Rosk ar et aHl201lh . care must be taken when computing the 
spiral phase of a pattern. One should select individual pat- 
terns by isolating a horizontal ridge (a single pattern speed) 
over a radial range where the signal significantly stands out 
from the noise, and proceed to calculate the spiral phase of 
that pattern. In each of the galaxies, we lift the most dom- 
inant patterns (the darkest horizontal ridges in the power 
spectra in the top panel of Fig. [5]) and look at the three 
quantities, VK c m (r,cj), jy s m (r,cj) and A m (r,uj) on the real 
and imaginary axis for each radial pixel in a ridge. We are 
able to calculate the spiral phase as: 

0?Jr,u>) = arctanf ^^^ ^ , (11) 

from which we can convert to real phase position within 
the domain to 2tt, by dividing by the multiplicity, m: 
p (r, Q) = 6™(r,uj)/m. Therefore, by considering only the 
Fourier coefficients of a single given uj whose amplitudes 
stand out from the noise in the power spectrum, we may 
retrieve 0JJ(r, uj) using equation JTTJ, an d calculate the cor- 
responding phase position of a single pattern as a function 
of radius and then the pitch angle. The phase positions for 
each case are shown in the bottom row of Fig. [5] The size of 
each point represents the amplitude, A m (r, cj), normalised 
to the largest amplitude of all the patterns shown in a given 
mode. Therefore the size of each point reflects the relative 
strength of a pattern at a given radius. The linear fitting 
of a constant pitch angle of equation © described at the 
beginning of the section is applied to points which exceed 
70% of the maximum amplitude for each pattern. We ap- 
plied a variety of these cuts and found that this cut best 
fit all the data. The best fit lines shown in Fig. [5] fit the 
points very well, and produce reliable pitch angles, whose 
averages are 12.9°, 20.3° and 33.1° for shear rates of ~ 1.3 



(simulation K), ~ 1.0 (simulation F) and ~ 0.7 (simulation 
R) respectively. This reproduces the expected trend of pitch 
angle vs. shear very well: the density wave mode patterns 
display more tightly wound structure when the shear rate is 
higher. 

3.2 Pitch angle of apparent spiral structures 

It is now interesting to look at the pitch angle of the appar- 
ent spiral arms in our simulations for two reasons: the ap- 
parent density enhancement of the spiral arms may manifest 
from more than one overlapping and interfering wave pat- 
tern (jGrand et al.ll2012al : IComparetta fe Quillenll2012h like 
the ones seen in the upper panels of Fig. [5] The other reason 
is that if the app arent spiral arms w in d with time (prev i- 
ously reported in iGrand et alJ [2012a.b: Wada et al. 201l|), 
the pitch angle must vary with time. This may have interest- 
ing consequences for observations of galaxy pitch angle vs. 
shear (Siegar et al. 2005, 2006), as winding spiral arms will 
show a range of pitch angles for a given shear rate, which 
will be reflected in observations of a galaxy sample. 

The method we use to trace the spiral arm peak position 
is a density weighting method, in which we select a point 
near the spiral arm of interest at some start radius (~ 5 
kpc), define an azimuth range that encapsulates the width 
of the spiral arm and weight by density to find the peak 
position (see [Grand et al.ll2012bl . for more details). This is 
iterated over a radial range until the spiral arm peak position 
is drawn out. This is done over a range of snapshots between 
1 and 2 Gyr of the simulation evolution and the pitch angles 
are calculated in the same way as the mode pattern pitch 
angles in Section 3.1. 

Because the pitch angles of spiral arms appear to de- 
crease as their evolution proceeds from birth to death i.e. 
the arm winds, care must be taken that spiral arms are not 
chosen at some special moment when they at a particular 
stage of winding. To ensure that we measure pitch angles in 
an unbiased way, only those spiral arms that show a single 
density peak in the azimuthal direction for each radius in the 
specified radial range for fitting are chosen for measurement. 
This is a limitation of this method. Fig. \6\ demonstrates an 
example of the application of the arm tracing criteria to one 
of the spiral arms in simulation K. Because it is possible 
to reliably trace spiral arms which show only single peak 
structure for the radial range considered for fitting, we re- 
ject those snapshots that show the spiral arm with messy or 
double peak structure, which typically occurs during spiral 
arm formation (t — 1.152 Gyr in Fig. [6]) and after the arm 
shows bifurcation or breaking (t = 1.2 Gyr in Fig. [6j). 

The results for several spiral arms in each simulation are 
shown in Fig. [71 It is obvious that every spiral arm pitch an- 
gle decreases with time. This is consistent with winding, co- 
rotati ng sp iral arms which have been reported in I Wada et al.l 
(|201lh and lGrand et al.1 (|2012al lbh. The mode pattern pitch 
angles calculated in the previous section are highlighted by 
the shaded regions in Fig. [71 The apparent pitch angle of 
any given spiral arm straddles the mode pattern pitch angle 
region calculated from the corresponding simulation. The 
overlap occurs close to the middle of the traceable time se- 
quence of a spiral arm. A snapshot of a time when both pitch 
angles overlap for each simulation is shown in Fig. [51 The 
agreement between the mode pattern pitch angles and the 
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t= 1.152(Gyr) 



t= 1.160(Gyr) 



t= 1.168(Gyr) 




200 300 





t= 1.184(Gyr) 



t= 1.192(Gyr) 



t= 1.200(Gyr) 






Figure 6. Snapshots of the disc density in polar coordinates. Density contours are overlaid in white. The traced spiral arm position is 
highlighted with a black line. Note the double peak structure at R ~ 5.5 and ~ 9 kpc at snapshots t = 1.152 and t = 1.2 Gyr, that 
prevents an unambiguous fitting to a single peak. This criteria defines the time range in which we can trace the spiral arm. 



apparent spiral arm pitch angles, together with the presence 
of several patterns of seemingly varying pattern speeds is 
consistent with the idea that two or more density waves may 
periodically constructively interfere with one another to cre- 
ate a spiral arm. The apparent pitch angle begins larger than 
the mode pitch angle, and approaches the mode pitch angle 
while the density grows (constructive interfering). The waves 
then pass and move away from one another, which leads to 
a stretch in the azimuthal direction of the apparent spiral 
arm density. This produces a winding in the overall apparent 
spiral arm pitch angle. The co-rotating spiral arm features 
seen in previous simulations may be understo od as the in- 
terferen ce of multiple wave modes as noted by I Grand et al.l 
(|2012aT ) who commented: "If there are indeed several wave 
modes present, it is evident that they must conspire in a 
specific way in order to produce a spiral arm feature that 
is apparently co-rotating" . It may be that the conspiracy is 
an interference mechanism of several wave mode patterns, 
each co-rotating at their respective radii of peak amplitude, 
that yie lds a near co-rotating pattern speed at the overlap 
regions (|Comparetta fe Gu illen 2012). The growth of such 
wave mode patterns is naturally expected from swing am- 
plification. However, the evol ution of such modes s hould be 
non-linear and complicated ([D'Onghia et al.ll2012T h which 
deserves further study. In this case, the wave mode patterns 
at the outer radii have lower pattern speeds than those at 
the inner radii i.e. AO/ AR < 0, so the evolution of the pitch 
angle may be represented by the time dependent gradient of 
the linear best fit line, 



b(t) = 



AR 



R(A0 o - AQt) ' 



(12) 



where AQ is the difference in pattern speed between different 
radii. Therefore, if AQ/AR < 0, b decreases with time, and 
hence the pitch angle decreases with time as well. 

The pitch angles calculated from the overall density en- 
hancement in Fig. [3 are plotted against galactic shear in 
Fig. [9] The range of the apparent pitch angle for a fixed shear 
rate appears to be similar to the scatter seen in the observed 
pitch angle-shear rate relation (e.g. Fig. 3 of ISeigar et al.1 
2006). This indicates that spiral arms are in fact really 
winding, possibly as a consequence of multiple mode pat- 
terns interfering, and that observers are seeing spiral arms 
at snapshots at varying stages of their evolution, and there- 
fore detect a range of pitch angles at a given rate of shear. 



4 CONCLUSIONS 

We have presented three N-body simulations each with dif- 
ferent shaped rotation curves, and analysed the pitch angle 
of spiral arms and underlying mode patterns with two dif- 
ferent techniques. We presented and compared the results of 
both techniques, and come to the following conclusions. 

(i) We calculate pitch angles of individual mode patterns 
for the first time by lifting their phase information directly 
out of the power spectrum of a strong mode multiplicity, 
and assuming logarithmic spiral arms in order to use the 
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t= 1.256(Gyr) 



t= 1.376(Gyr) 



t= 1.176(Gyr) 




Figure 8. Face on view of each simulation (from left to right: simulations R, F and K) when the apparent spiral arm pitch angle coincides 
with the calculated mode pattern pitch angle. The degrees of spiral morphology are strikingly different: spirals become increasingly tight 
going from left to right. 
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Figure 7. Pitch angle evolution of the apparent spiral arm for 
simulation R (circles), F (triangles) and K (diamonds). In all cases 
the pitch angle decreases with time, which indicates the winding 
nature of the overall density peak. The shaded horizontal regions 
represent the mode pattern pitch angles, determined from the 
patterns in Fig.[5]for simulations R (red), F (green) and K (blue). 
Note that the range of apparent spiral arm pitch angles clearly 
map out separate domains and straddle the mode pattern pitch 
angles of their respective galaxies. 



linear fit to get a constant pitch angle from the gradient. 
We find the mode pitch angle is correlated with the rate of 
galactic shear: the pitch angle is smaller for higher galactic 
shear rate and vice versa. This is consistent with the an- 
alytical predictions bas ed on swing ampl i ficatio n theory in 
ljulian fe Toomrd (| 19661 ) and lFuchs et~al1 (j2005l ). 

(ii) Pitch angles of the spiral arm density enhancement it- 
self are calculated in the same way from the density weighted 
position of the spiral arm peak. These pitch angles decrease 
with time, as the spiral arms grow from a relatively open 
morphology, then wind over time to become more tightly 




Figure 9. All calculated apparent spiral arm pitch angles plotted 
as a function of galactic shear for simulations R (red), F (green) 
and K (blue). 



wound until they disrupt. This is consistent with previ- 
ous si mulations that reported winding and co -rotating spiral 
arms (|Wada et al.1 1 20 111 : I Grand et alj|2012el lbh. 

(iii) The range of the variation of the pitch angles owing 
to the winding spiral arm features are similar to the scatter 
seen from the observed relation between the pitch angle and 
the shear rate in spiral galaxies. This indicates that real 
galaxies exhibit transient, winding spiral arms that may be 
caused by interfering spiral waves. 

The pitch angle is a key structural parameter of spiral 
arms, which is shown to depend on the shape of rotation 
curve. In the context of both density wave mode patterns 
and the overall stellar density enhancement, the morphology 
of spiral arms is important to the evolution of disc stars that 
themselves give rise to the structure. Phenomena such as 
radial migration is likely to be dependent on the orientation 
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of the spiral arm against the orbit of stars, and hence the 
spiral morphology of the galaxy in general. The pitch angle 
of the spiral arm may therefore affect how much a star can 
migrate. It is possible that there are other criteria which a 
star must satisfy that, coupled with the spiral morphology, 
will determine the ability of a star to migrate radially. These 
criteria will be explored in a forthcoming paper. 
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